Growth of perturbations in dark matter coupled with quintessence 



Tomi KoivistcQ 

Helsinki Institute of Physics, FIN-00014 Helsinki, Finland and 
Department of Physics, University of Oslo, N-0316 Oslo, Norway 
(Dated: February 2, 2008) 

We consider the evolution of linear perturbations in models with a nonminimal coupling between 
dark matter and scalar field dark energy. Growth of matter inhomogeneities in two examples of such 
models proposed in the literature are investigated in detail. Both of these models are based on a 
low-energy limit of effective string theory action, and have been previously shown to naturally lead 
to late acceleration of the universe. However, we find that these models can be ruled out by taking 
properly into account the impact of the scalar field coupling on the formation of structure in the 
dark matter density. In particular, when the transition to acceleration in these models begins, the 
interaction with dark energy enchances the small scale clustering in dark matter much too strongly. 
We discuss also the role of an effective small scale sound speed squared and the issue of adiabatic 
initial conditions in models with a coupled dark sector. 
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I. INTRODUCTION 



Cosmological observations suggest that about ninetysix percent of the present energy density in the universe resides 
in forms unknown to the standard model of particle physics. Measurements of the fluctuations in the cosmic microwave 
background (CMB) radiation are consistent with a nearly critical density universe^ , as predicted by the early inflation. 
Assuming the validity of general relativity, roughly two thirds of this density nowadays must be due to a negative- 
ly — pressure component called dark energyJ2| (DE) to explain the acceleration of the universe indicated by distance-redshift 
. curves of the high rcdshift supernovaej3|- The rest of the density is in gravitationally clustered matter with negligible 
pressure on cosmological scales. However, only a fraction of this matter consists of the (partly luminous) baryonic 
matter known to the standard model. Dark matter (DM) is needed for consistency with a convincing number of 

■ astrophysical observations, for example the flat rotation curves of galaxiesQ. 

It has been proposed that the two dark components do not evolve independently but are non-minimally coupled. 

■ Due to their different pressures, the densities of these components decay with different rates. This rises the disturbing 
Or question why they appear to be of the same order of magnitude right now. An explicit coupling between DE and 

DM 5] could provide a phenomenological solution to this so called "cosmic coincidence problem" 0, 0, H, EJ- On 
the other hand, there is theoretical motivational to couple quintessence fields^], ^2 to DM (see, however |13|). In 
particular, if it is held that scalar-tensor theories should be described in the Einstein frame, it follows that scalar fields 
with a nonminimal coupling to other matter are general predictions of all dimensionally reduced theories |l4l[T5j . In 
, the context of string theory, the dilaton could be such a field, and indeed has been proposed as a candidate for DE 
0,0]. Since generally the dilaton coupling to different matter species could be different [T^ and such couplings to 
baryonic matter arc tightly constrained by high precision solar system tests of general relativity 19], it is of interest 
to concentrate on simpler models where the dilaton is only coupled to DM, while the possible coupling to baryons is 
assumed negligible. 

This paper concerns the linear perturbations in such models. The purpose is to find how the formation of cosmo- 
logical structure in pressureless, dust-like matter is influenced by a coupling to a scalar field, when this coupling is 
significant already at the background order as required if it is to address the coincidence problem. First the cosmo- 
logical field equations for scalar field interacting with DM are derived to linear order in section II. These equations 
have been presented in the literature before (for example, see |20j). and the reader familiar with them may skip 
this section. However, there we hope also to somewhat clarify the relations of different possible formulations and 
theoretical frameworks for the coupled dark sector. In the two following sections we investigate specific cases. The 
simplest models of interacting DM and DE can be excluded by comparison to the CMB observations [2lJ, but more 
sophisticated scenarios exist 0, 0, 0|. In section III we study a model featuring two DM components, only the 
subdominant of them interacting with quintessence |22j. A model where the DE-DM coupling is given by a polynomial 
function of the DE field0] is considered in section IV. These two models are based on a low-energy limit of effective 
string theory action and have been previously shown to naturally lead to a late-time acceleration of the universe. 
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Though the details of the mechanisms are different, the DE domination in each case is triggered by the time-evolving 
coupling. Large scale structure in both of the two scenarios have previously been claimed to fit observations as well 
as or better than the standard ACDM model. However, here a closer inspection shows that these models are ruled 
out because of the excessive DM structure they predict in the late-time universe. In subsection IVB we consider the 
adiabatic initial conditions in for a scalar field coupled to DM. In section V we discuss the relation of our results to 
other models and conclude. 



II. EQUATIONS 



A. Field equations and Bianchi identities 



We consider an action of the form 



(1) 



where \ stands for an interacting matter species, and m for non- interacting ones. The constant K = (8ttG) 1 ^ 2 . We 
use the definition 



(2) 



and the Lagrangian of the field cf> is taken to consist of the two first terms in the action integral Eq. • The Einstein 
equations resulting from this action can be written as 



and the equation of motion for the field <p is 



The Bianchi identities then imply 



□0 = f'(0)-7^)/; (x) . 
v^ i(m) = o, 



(3) 
(4) 
(5) 



V^to = (J^Cx) _ T M(x))2xg) (6) 

7xW 

Thus the energy momentum tensor for the interacting components are not separately conserved. 
From now on we will assume that 

£M=-Vm4T-4)) ^- , (7) 

where \x l \ <C 1, i.e. the interacting component behaves as dust. This is in fact the case most commonly studied in 
cosmology (for models with DE coupled to neutrinos, see Refs. [23ll2^| '). One notes that the term j x ((f>)£^ in the action 
Eq.Q can be naturally interpreted as a Lagrangian for particles with field-dependent masses, m a (<p) = m a 7 x (</>). Such 
a Lagrangian can follow from a conformal transformation. Consider a scalar-tensor theory defined by a gravitational 
Lagrangian proportional to R/ \f~f{<t>)- Then the metric in the Einstein frame is <7 M „ = g^v/ \/l{<j))- In terms of this 
metric, the action for the Lagrangian Eq.Q reads 

s x = J tfxjqHwtM . (8) 

(At this point one may have to redefine the field in order to have the canonical kinetic term for the scalar field as 
in the action Eq.Q.) A more standard way, however, to describe scalar-tensor theories in the Einstein frame is to 
define the component 

rW=7x(0)TW. (9) 
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This energy momentum tensor then depends on the field (j>, which turns out to be coupled to the trace T^ x \ For some 
purposes the ^-formulation is more convenient, since there the divergence in Eq.© vanishes except for the spatial 
part in the linear order. However, we will also use the X-formulation for some calculations, especially in subsection 
IIIB. With the variable mass particle interpretation, the difference is simply that p x is proportional to the number 
density of the DM matter particles, while px is their gravitating energy density. The fact that the action Eq.0J 
is equivalent to the case of a trace-coupled scalar field and the possibility to make direct contact with conformally 
transformed scalar tensor and variable mass particle theories are independent of the form of the coupling. However, 
it is required that the Lagrangian is effectively of the form Eq. JJJ . From that follows (when | x % | <C 1 ) 

£(x) = T (x) = _ px _ ( 10 ) 
These simple identities would not, of course, hold for a more general Lagrangian. 



B. Cosmological background and perturbations 

We consider a flat FRW universe. Including small perturbations about the background, we can write the line 
element as 

ds 2 = a 2 (T)[-dr 2 + (S i:j + h tJ )dx % dx j ]. (11) 

This parameterisation of the metric perturbations is called the synchronous gauge. There we can characterize the 
scalar modes of hij in the Fourier space as 



kikjh(k.,T) + (hkj - ~Sij)6r](k, r) 



(12) 



when k = kk. In the following, prime is a derivative with respect to <p, an overdot represents derivative with respect 
to conformal time r, and H is the conformal Hubble parameter H = a/ a. These and mostly all other notations and 
conventions are those of Ma and Bertschinger 25] . With these, the Friedmann equation reads 

H2 = &i(& + VW+1 x p x +Pm )' (13) 

while evolution of the scalar field is governed by the equation 

4> + 2Hj ) + a 2 V(0)' = -a 2 y x (0) Px , (14) 

following from Eq.Q- In addition to the interacting species and the scalar field, we include standard CDM, baryons 
and radiation in the universe. Then the energy densities scale as usually, and therefore, with the scale factor normalized 
to unity today, we may write 

0° n° n° n° 

Pc Pb Pr firs 

PX = "J > Pc = -3, Pb = -3, Pr = —■ (15) 

a J a A a A a 4 

Here superscript denotes the density today. However, instead of p x the gravitating energy density is now the density 
px = IxPx fro m definition Eq.Q, as implied by Ea.l |13(l . 
The equations for the perturbed metric are 

k 2 V -±Hh=^6T°, (16) 
k 2 r) = ^- 2 ih8Tl (17) 



h + 2Hh-2k 2 r] = ^7?, 



(18) 
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h + 6ij + 2H(ti + 677) - 2k 2 v = —(hkj - -<%)£j, (19) 



where the components of the perturbed stress energy tensor are 

ST° = -p r S r - Pxil'x^ + lx$x) _ Pb$b - pc$c ~ ~^H> ~ Vf, (20) 

ik t 6T° = ^p r S r 9 r + l x p x 9 x + p b b + <ftk 2 ip, (21) 

Stf=pr5 r +34pbSb + ^4xp-3V'<p. (22) 

We have grouped together photons and neutrinos with subscript r, since their contribution is as usual, and used 
(p to denote 6cj>/<f) as in 26] . In addition to the field equations, we need the equations of motion for the individual 
components. For the perturbations of the coupled scalar field, we derive the Klein-Gordon equation from Eq.Q to 
find 

(p + 2Hp + {a 2 V" + k 2 )ip + X -4>h = -a 2 p x (7^x + 7»- (23) 
For the x-component, the Bianchi identities Eq.© in the linear order are 

s x = -e x - h, (24) 
e x = -Hd x + ^{k 2 tp - #y . (25) 

7x 

Conservation equations for uncoupled matter stay unmodified from the ones given inplj. 

III. MODEL WITH TWO COLD DARK MATTER SPECIES 
A. Background solution 

In this section we investigate the model proposed by Huey and Wandelt|22j|. In this model quintessence with an 
exponential potential, 

V{ct>) = V e' Ka t (26) 

is exponentially coupled to x- m atte r 

7x (0)=e re W-«. (27) 

There is also present an uncoupled species of cold DM. At early times, the field behaves as standard exponential 
quintessence, but when it reaches the minimum of its effective potential, V e ff — V + ^ x p x , it evolves to a potential- 
dominated stage driving the acceleration of the universe. 

There is motivation for the additional complexity of introducing two kinds of DM. It is well known that with 
minimally coupled exponential quintessence the tracking solutions do not yield an accelerating phase p^|. If the field 
is exponentially coupled to matter, there exists an accelerating attractor|23, but obviously the field could not have 
reached it until today lest there would not have been a matter dominated era at all. Moreover, it is known that 
in the such models the coupling between DE and DM is tightly constrained by the CMB spectrum|2!| and other 
cosmological data[3(| . 

The amount of quintessence during BBN must be small, constraining a > 5 |27| ■ The onset of acceleration should 
occur after matter-radiation equality, which is achieved when n<p c > 108/a. Finally, consistent late-time background 
evolution requires a/ (3 ^ 3/5 |22j . As a representative parameter combination we will use a = 5, (3 — 15, 4> c = 25 in 
numerical calculations. Fig.lJTJ shows then the evolution of the fractional energy densities f2j = k 2 pi / (3H 2 ) . We have 
arranged the interacting and standard CDM components to be of equal importance today. We see some transient 
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oscillations at the late epoch when the coupling term begins to hold the field at the minimum of its effective potential. 
However, asymptotically the coupling will keep the ratio of quintessence and interacting DM densities constant and 
the field evolution is again characterized with a constant scaling. 

Consider the component px = IxPx- Since during the tracking phase, assuming a constant background equation 
of state wb, 



3(1 + wb) , 
K(p = log(a) + constant, 



(28) 



then px ~ a 3 ' 1 +'q( 1 + 10b )1. Thus this energy density is growing in a phantom-like manner. Its effective equation-of 
state parameter, defined by 



px + 3Hp x (l + w^ ff) ) = 0, 



(29) 



is shown together with the total equation of state of the universe in the left panel of Fig.JSJ). Similarly we define the 
effective equation of state for the field <j>: 



p + 3ffp (l + ^ e//) ) = 



(30) 



Then we have from Ea. lfH|) that 



P4, + ZHip^+p^) = -i {<t>)p x 4> 



w 



(eff) _ P<P+l X PxH(^ H ) 



(31) 



We plot both u>0 and in the right panel of Fig.®. 




FIG. 1: Evolution of the fractional energy densities in model of section III. The solid line is fi^,, the dashed line Qx, the dotted 
line is fi r and the dash-dotted line (Q c + Qt)- 
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FIG. 2: Left panel: Total equation of state w (solid line) and the effective (dashed line) for interacting CDM in the 

model of section III. Since w < —1/3 today, the universe accelerates. Right panel: Equation of state for the quintessence in 
the same model. Solid line is the formal w^, and dashed line is the effective w^^\ 



B. Perturbation evolution 

Defining C = log(7 x ), the coupled perturbation equations for the component X can be written as [2? 



ox 



9 x --h + Ctp + C"<lxp, 



(32) 



?x = -H9 x + c'(k 2 <p-4>e x ). 



(33) 



For the largest scales the source term in Ea. (|33|l can be neglected, and thus no velocity perturbations are generated. 
Since then S c — Sx — C'ip, the coupling has a small effect on Sx, which then follows closely the standard CDM 
overdensities. 

To find the small scale evolution of perturbations, we follow derivations found in Refs. [2fll2^ . l3lll32] |. Using a linear 
combination of Eas. (|16fl and (|18fl . 



h + Hh=^ [ST° - STf] , 



(34) 



and averaging over the oscillating time derivatives (p, ip and S r , one can combine Eqs. (|32() and 1|33|) into 



S + HS 



2k< 



PxSx - C'W = i-C'k 2 + HC"<\> + C - —V'] V + ^{pcSc + p b Sb), 



when cl ~ 0. Averaging out the small and oscillatory (p and (p also from Ea. (|23|l yields 



a 2 p x C'S x + \4>h 
a 2 (V" + Px C") + k 2 



(35) 



(36) 



This approximation is in fact reasonably good for all scales. We have tested it by comparing the LHS and RHS of 
this expression evaluated by numerically integrating the full perturbation equations^- Inserting the approximation 
in Ea. H35|) and using Eq. (|33|l . we get in the small scale limit 



rt\ : // < )S X = \U 2 



C' 

n x ( i + 2— ) 5 X + n c Sc + ftbSb 



(37) 



Here we have used C" = 0. During tracking era, Cl x ~ a 3 ^ 1+WB ^a + 1 ) '•I vanishes to a good approximation. If we 
assume that then the noninteracting matter perturbation scales as S c ~ r m ~ 2 , we can use the tracker solution Ea. l28|) 
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to find that 



Sx + 



33 , 

1 + — {l + w B ) 
a 



-S x = Ar m , (38) 



(1 + 3w b )t 

with the constant A determined by the magnitude of noninteracting matter perturbations at a given time. The full 
solution to this equation is 

5 x =D + B T 1 - b + - 4t? 7^r m+2 , (39) 

(b + m+l)(m + 2) y ' 

where B and D are integration constants and b — [2 + 6/3(1 + wb) /a]/(l + 3wb) ^S> 1. Thus the last term gives the 
dominant mode, which has the same scaling as the noninteracting DM. Therefore the ratio Sx/S c is constant during 
tracking era. During matter domination, when m = 2, this means that 5x ~ a, and Sx/5 C = 3(1 — 3/a 2 )/(10 + 
128/a) < 1. 

When the field reaches the minimum of its effective potential, it evolves with the ratej2^ 

4 = ~^-sH. (40) 
a + p 

The critical period we are interested in is the transition to the joint domination of quintessence and interacting 
matter that is ongoing today. The field derivative is then oscillating about the value given by Ea. (|40|l . to which it 
will asymptotically settle (see Fig.©). Since Qx is now approaching Q c , and (1 + C' 2 ) ^> 1, we approximate Ea. ll3*7|) 
by dropping the source terms due to noninteracting matter. This approximation is well justified since, as it turns 
out, the RHS of Ea. H37l) is soon dominated by the rapidly growing perturbation in the density of interacting DM. We 
therefore write 

5x+{l + ^-^jH5x = \H 2 n x {l + 2p 2 )5 x . (41) 

To solve this equation analytically, we would need expressions for the background quantities H (which follows from 
Wb) and £lx- To that end we use the asymptotic value also for the latter 

w B =-0/(a + 0), n x = a/(a + 0). (42) 

Inserting these in Eq.|@IJ yields the solutions 

1 (a + 10/3) 



n ■■ 



4 (a + 0) 



-1± Wl 



24a(a + /3)(l + 2(3 2 ) 
(a + 10/3) 2 



(43) 



The exponent of the growing solution (corresponding to the positive sign) is large for typical values in this model. 
For a = 5 and 8 = 15, n w 11.2. This is in strong contrast to standard DM perturbations growth rate of which is 
slowing after the matter domination when 5 C ~ a. Therefore 

Ox-tc-Sx-aT, r = n-^-^- (44) 

2(a + p) 

which is, consistently, the growing solution of Ea, 133|l in this limit. Numerical results are shown in Fig.® 48] . 
There we see that the interacting DM density contrast adapts to this asymptotic scaling in fact as the transition to 
acceleration begins, although the solution Eas. (|40|l . (|42|l has not been reached. 

We conclude that when the energy in the x _m atter is becoming a significant contributor to the background expan- 
sion, the perturbation growth in this component is almost explosively sped up by the interaction pumping energy 
into it from the 0- field. As seen from Ea. l|3()l) . the DE is then inhomogenenous on small scales. It has been found 
also that in the nonlinear regime DE tends to clump with DM if coupled[H The large coupling for which 
C = (3 in this model results in a large effective gravitational constant in Ea. l41f> . G e ^ = (1 + 2/3 2 )G. In earlier 
cosmology this is not as important, since then the standard DM overdensities constitute the dominant gravitational 
sources. There the modified Hubble friction term in Ea. H37(l acts to reduce the relative magnitude of Sx to S c , but 
these evolve with the same rate. However, in the late universe the cosmic structure is drastically increased, except at 
the largest scales where we see only little oscillations. Consequently, the shape of the matter power spectrum in this 
model is highly incompatible with observations. Also, as the increase in structure will cause rapid deepening of the 
gravitational wells, the resulting CMB spectrum is modified by a strong late integrated Sachs- Wolfe effect due to time 
variation of the gravitational potentials. We have confirmed these considerations by numerical integration for several 
parameter values satisfying the relations given in the previous subsection. As a non-negligible amount of ^-matter 
and a relatively strong coupling are required in this model, we feel it safe to state although having not performed an 
extensive exploration in the space of all cosmological parameters, that this model cannot be succesfully confronted 
with observations. 
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FIG. 3: Interacting DM overdensities in the model of section III. The slope s of the curves in the logarithmic plot is s » 2, 
1, 11.2 early in the radiation domination (when the scale is outside the horizon), during matter domination and during the 
transition to accelerating universe, respectively. The fc-values from bottom to top are k = 2.4 • 1(T 6 , 2.4 • 10" 5 , 2.4 • 1(T 4 , 
0.0024, 0.024 and 0.24 Mpc" 1 . 



IV. MODEL WITH POLYNOMIAL COUPLING 
A. Background solution 

The model proposed by Bean and Magueijofl^ is another example of a modified exponential quintessence scenario. 
Again a field is coupled to DM in such a way that quintessence is dominant today and driving the acceleration of the 
universe while the model still has its parameters of the order of the Planck scale and its outcome widely independent 
of the initial conditions for the field. This is achieved with only one DM species by using a non-exponential form for 
the coupling. To be definite, we have in the action Eq.JQ) 

V(4)=V e-* a +, 7x(</>) = l + AK0-0o)] b . (45) 

For numerical values we will use a — 8, b = 8, A = 50 and K(po — 32.5. The background evolution is shown is Fig.|@J. 
In this model the field leaves the standard exponential track long before matter domination. After initial transients, 
the field tracks the background radiation, until the coupling term pushes the field into kination. This lasts until the 
field reaches the minimum of its effective potential. Then it will at due time drive the acceleration. However, as 
the matter density eventually dilutes, the coupling force will fade, leading the field back to standard scaling in the 
asymptotic future. Thus matter domination is restored and the acceleration in this model is transient. We plot the 
effective equation of state parameter for the field cj), as defined in Ea. i|31[l . in the left panel of Fig.©. In the right 

panel we plot the total equation of state for the universe together with the as defined in Eq.J55J). Since in 

this model px governs the background expansion also in the matter dominated era, it is important that the coupling 
term has then negligible effect on the scaling of this density, since a significant deviation from a ~ r 2 would result in 
modified conformal distance to the last scattering surface, possibly dangerous to the compatibility with observations 
of the peaks in the fluctuation spectrum of the CMB radiation j^. One also expects that then the growth rate of DM 
perturbations will not be much modified. However, one would guess that the nonstandard evolution during transition 
to DE domination, seen in Figs.(0J and (J5J, has some effects on the perturbation growth as was the case in the model 
of section III. Indeed these effects prove fatal to this model, as will be seen below. 
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FIG. 5: Left panel: Total equation of state w (solid line) and the effective w^' (dashed line) for interacting DM in the model 
of section IV. Since w < —1/3 today, the universe accelerates. Right panel: Equation of state for the quintessence for the same 
model. Solid line is the formal m^, and dashed line is the effective wi e . 



B. Perturbation evolution 



The linearized Einstein equations have in fact been solved in this model 36]. However, this was not done taking fully 
into account the impact of the coupling to the matter perturbations, but instead by assuming that the fluctuations 
in p x decouple from those in (f>. This was realised by setting the DM velocity perturbation to zero, as conventional in 
the synchronous gauge. However, now the interaction with DE will drive 9 X to evolve even if its initial value is chosen 
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to be zero. One sees that the source term in the Euler equation Ea. (|25|l vanishes only in the absence of coupling or 
if the DE is perfectly smooth. Neither is the case now. We show the effect of physical assumptions in Fig.©. 

Since now C" = 0, Eq.(j33 is not valid. However, using the approximations that lead to that equation, we can 
derive the corresponding evolution equation for the quantity S x . It should capture the qualitative behaviour of the 
DM perturbation, since Sx = S x + C'tp « S x for small scales. The equation we find reads 



H 



6-1 



1 + Xy 



-K(j) S 



-H 2 n x 



(1 + Ay") 2 



(46) 



where y = «(</> — <fio)- The effective gravitational constant can be thought to be the expression is square brackets 
times G. We plot G^-^ /G — 1 as a function of 4> in Fig.Q for a choice of parameters. The values of G e ^ becomes 
relevant at the matter dominated era, when the source term in the RHS of Ea. (|46() is non-ncgligiblc. However, as 
the coupling has driven (f> near to (f>o, we have then G^L ~ G (for our example parameters, —0.5 < y < from 
a = 10~ 4 until a — 0.1). This changes when fl^ climbs up, since then tfi rolls forward. Similarly as in the previous 
case, Sx is driven to rapid increase. We show its evolution in Fig.JHJ, where we also plot the baryon overdensities 5b 
for comparison. The boost experienced by Sx is not as vast as we saw in the previous section, but still results in an 
order of magnitude or two too large Sx- Again we believe that notable improvement is not found by adjustments of 
parameters, but by revision of the scenario. 




FIG. 6: The CMB spectrum in the model of section IV (solid line). Dash-dotted line corresponds to suppressing the coupling 
of dark energy and DM perturbations (i.e. 9 X = 0), and dashed line corresponds to suppressing DE perturbations (i.e. ip = 0). 



Finally, a remark about the adiabatic initial conditions we have used in the calculations. Adiabaticity of 
perturbations [^3] means that given any two super horizon regions A and B in a perturbed universe at an instant 
r, B is identical to A at some other instant r + St. Evolution of the regions is the same, but slightly asynchronous. 
Then the ratio of fluctuation and time derivative of say a density field is common to all density fields. Thus we can 
define the entropy between species i and j as 
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FIG. 7: The deviation of the effective gravitational constant from the standard one, Gr a ' * ' /G — 1, as a function of the field K(f> 
for the model of section IV. The solid curve corresponds to the Gjjf m Eg . 1461 . and the dashed curve is 2C' 2 , corresponding 
to G e J^ as it appears in Eg. 13711 . For both definitions the deviation is small when <f> m <f>o, but becomes very rapidly large when 
the field <f) evolves to larger values. 
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FIG. 8: DM overdensities (solid curves) in the model of section IV. Dashed curves are Sb- The fc-values from bottom to top 
are k = 2.4 • 10" 6 , 2.4 ■ 10" 5 , 2.4 ■ 10" 4 , 0.0024, 0.024 and 0.24 Mpc" 1 . 



The general adiabatic condition then dictates that all Sij 's together with their derivatives vanish. Note that only in 
the absence of interactions does this condition coincide with the usually cited 



1 + Wi 1 + Wj 



(48) 



There can be difference between these conditions, as seen for the cases we are now considering in Figs.(j2J and JSJ. 
For a scale far outside the horizon and deep in the radiation dominated era, the growing adiabatic mode |25l |38j of 
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the scalar field perturbation satisfies 

30 2 - (iy - h' xPX )4>r + \ Px [ 7 > 2 - iyv - i 2 a 2 Px 



(6(f> + j' a 2 p x T)V - Y'p x <j) 2 T 



-<j)Ak 2 T 2 , (49) 



(60 + 7^a 2 p x r)F' - i{p x ^ 2 r 

Normalization A = —1/2 is convenient and sets comoving curvature perturbation equal to unity. Fortunately the 
above expressions can in many occasions be simplified by dropping the square bracket terms, but we have kept them 
here for the sake of completeness. The departure of the usual adiabatic condition Ea. (|48|l from the generalized one 
has been pointed out also in Refs. [39L ficj . If the field is a tracker, the initial conditions may indeed be irrelevant also 
for the perturbed field. However, for minimally coupled fields it has been shown that in some instances this is not the 
case 0,^2) while isocurvature perturbations in coupled quintessence remain to be systematically studied. 

V. DISCUSSION 

Interestingly, in the two models we have investigated, the growth of DM perturbations is enchanced, whereas in 
minimally coupled models the late time effect of DE is the opposite. This may be understood by comparing to the 
case where the interacting DE and DM satisfy the adiabatic relation S^ x = S$ x = 0, where S,p x is defined by Eq.l|47|l. 
Then the system of these components can be treated as a unified fluid without intrinsic entropy. It is known that 
the effective sound speed squared c| = p/p of such a fluid governs the late dynamics of the perturbations 0|. If 
c| > 0, they are driven to oscillate, if c| < they tend to blow up. This is the work of the pressure gradient sourcing 
the evolution equation one can derive for the unified fluid density contrast, 

S+[l- 3(2w - c 2 s )]H5 - ^(1 - 6c| + 8w - 3w 2 )H 2 S = -k 2 j- p 5 - ( 51 ) 

This equation is exact (to linear order, in the comoving gauge). As w and c| are given by the background expansion, 
one needs only information of the pressure perturbation to solve the structure in the universe. In the adiabatic case, 
Sp/Sp — Cg, but in the scalar field models considered here this assumption is not valid. However, using equations 
from subsection IIB, the small scale limit of the pressure gradient in the late universe is readily found. It is 

- k T P = - k p-^wm m (<Wi ( ] 

(also in the comoving gauge) and turns out positive for the models of section III and IV. Thus one may regard the 
enchanced small scale growth rate there as a consequence of negative effective (non-adiabatic) sound speed squared. 
Now the pressure gradient approaches a /c-independent expression and the growth rate saturates at a given value when 
k is increased beyond some critical scale, unlike in the adiabatic case. Note also that the total density perturbation 
5 in Eq.JSTJ, 

s = Pxi' x w$x + + vwiv + g# (53) 

is directly relevant to the integrated Sachs- Wolfe effect in the CMB fluctuation spectrum. 

Before concluding, we use the opportunity to comment on yet another model with DE-DM interaction 0- There 
DE with a constant equation of state parameter is coupled to DM in such a way that 

p^ + iHp^ = 3HX 2 ( Px +p^), (54) 

px + 3H Px = -3HX 2 { Px + pt), (55) 

where A 2 is a constant measuring the intensity of the interaction. It might be difficult to find a field theory underlying 
this phenomenological model, but that would be necessary to be able to derive the corresponding equations to the 
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linear order. Otherwise they must be put in by hand. The linear DE continuity equations presented in the recent 
paper^5| imply that there 

= 2,H P4> \ 2 {8^ + r5 x ), (56) 

V„Tf W = [3H P4 ,\ 2 (l + r)0 ],i, (57) 

where r = p x /p r f > . Except perhaps for the sign of the 0-component of the divergence^?!, these expressions are not 
unreasonable, but do not follow unambiguously from Eas. l|54f) and (|55|l . However, once the divergence of the DE 
stress tensor is given, the interacting DM continuity equations are uniquely fixed. When Eas. l|56|l and 157|) hold for 
the DE and the DM Lagrangian is effectively of the form (JJJ) , the Bianchi identities determine that 

Sx = -e x -\h-2,H\ 2 {h <t> + 8 x ), (58) 



9 X = -H9 X + 3ffA 2 (l + -)6a. (59) 

r 

From these we recover the equations used in Ref.^]| only by setting 9 — in the first equation and forgetting about 
the second one. But as detailed in the previous section, such approximation can leave crucial features out from the 
evolution of perturbations. Therefore we must regard that the results of the data analysis in Ref.^j|, where it is 
announced that observations prefer interacting DE models, call for reinvestigation. 

As mentioned in the beginning of this paper, the origin of most of the energy density in the universe is unknown. 
Suggestions can be based on speculations of string theory, but this is tentative at the present. However, some inventive 
mechanism, as the ones considered here, is necessary if some light is required to be shed on the coincidence problem. 
Moreover, as the results of this paper demonstrate, considerations of cosmological structure provide powerful means 
to constrain such mechanisms. 
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